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We discuss the possibility that the X(3872) can be a hadronic DD* bound state in the framework 
of a potential model. The potential is generated by the exchange of pseudoscalar, scalar and vector 
mesons resulting from the Lagrangian of heavy hadron chiral perturbation theory. The hadronic 
0\ ' bound state configuration contains charged and neutral DD* components, while orbital S- and D- 

waves are included. Isospin symmetry breaking effects are fully taken into account. We show that 
binding in the DD* system with J PC = 1 ++ already exists for a reasonable value of the meson- 
exchange regularization parameter of A ~ 1.2 GeV. We also explore the possibility of hadronic BB* 
bound states and show that binding can be achieved in the isoscalar limit for J PC = 1 ++ or l" 1 . 
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I. INTRODUCTION 

The new resonance X(3872) has been discovered a few years ago by Belle [1] and has been later confirmed by the 
CDF [2], DO [3] and BABAR [4] collaborations. It is classified as an isosinglet state with positive charge parity. Current 
averaged results for the A(3872) mass and width are: m x = 3872.2±0.8 MeV and T x = 3.0±^±0.9 MeV [5]. For the 
i—i, A(3872) several structure interpretations have been proposed in the literature (for a status report see e.g. Refs. [6-8]). 
In the context of molecular approaches [9]-[33] the A(3872) can be identified with a weakly-bound hadronic molecule 
\ whose constituents are D and D* mesons. The reason for this natural interpretation is that rax is very close to the 
D°D* threshold and hence is in analogy to the deuteron — a weakly-bound state of proton and neutron. Note, that 
the idea to treat charmonium-like states as hadronic molecules traces back to Refs. [9, 10]. Originally it was proposed 
that the state A(3872) is a superposition of D°D*° and D°D*° pairs. Later (see e.g. discussion in Refs. [17, 19, 20]) 
also other structures, such as a charmonium state or even other meson pair configurations, were discussed in addition 
to the D°D*°+ charge conjugate (c.c.) component. The possibility of two nearly degenerated A(3872) states with 
positive and negative charge parity has been discussed in Refs. [26, 34]. 

In several papers [12, 29, 31, 33] the possibility that the A(3872) can be interpreted as a hadronic molecule — a 
■ bound state of D° and D*° mesons — has been investigated in potential models. These approaches can be traced 
back to Ref. [11], where possible deuteron-like two-meson bound states have been considered in the context of a 
potential generated by the pion exchange mechanism. It was shown that bound states of two mesons are possible in 
analogy to a deuteron model with one-pion exchange only. In addition, the relevance of the tensor interaction has been 
pointed out [11] which connects different spin-orbit configurations. These tensorial terms were found to be important 
to generate bound states in the potential model of the deuteron with a reasonable value for the cutoff A, which is 
a parameter regularizing the hadronic interaction. The value of A is determined phenomenologically and depends 
on the model, though its scale in low-energy hadron physics is of the order of 1 GeV. Later, after the discovery of 
the A(3872), Tornqvist also showed [12] the importance of isospin breaking effects in a possible binding mechanism 
because of the mass differences between the neutral and charged D(D*) mesons. Thomas and Close reinvestigated [31] 
in the context of one-pion exchange if the A(3872) can be a J p = 1 + DD* bound state including charged mode and 
D-w&ve configurations. They also point out technical differences with respect to previous approaches, but moderate 
binding depends rather sensitively on the parameters involved. However, Liu et al. [29] claim that when taking into 
account both pion and also sigma meson exchange potentials the D and D* mesons cannot form a bound state with 
reasonable values for the potential parameters. They show that a bound state is not present for values of A < 5.8 
GeV, where the upper value is too large compared to the typical hadronic scale of A ~ 1 GeV. But in their calculation 
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only the S'-wave configuration was taken into account. Moreover, only the interaction transitions D°D*° — ► D*°D° 
involving neutral D mesons were considered. This causes that the interaction strength is reduced by a factor of 1/3 
compared to that of the isoscalar state in the isospin symmetry limit. Later, Liu et al. [33] showed that the A(3872) 
is open to the possibility of a loose bound state of DD* when one considers the isospin symmetry limit with 1 = and 
also when further heavy meson-exchange, such as p or u exchange, are included. Again, they restricted the calculation 
to the S-wave component and to the isospin symmetry limit. 

The main objective of the present work is to improve the description of the A (3872) as a possible bound state 
of D and D* mesons in the context of the potential model considered previously in Refs. [11, 12, 29, 31, 33]. In 
particular, we first consider a full set of mesons (it, a, rj, p and to) in constructing the nonrelativisitic one-meson 
exchange potential, then we include charged and neutral DD* components both in orbital S- and D-waves. Full 
account is taken of isospin breaking effects by incorporating mass differences in the D and D* mesons. As in the 
preceding papers we also regularize the potential by introducing form factors containing the cutoff parameter A. 

In the manuscript we proceed as follows. First, in Sec. II we discuss the basic notions of our approach. We define the 
effective meson Lagrangian based on heavy hadron chiral perturbation theory, which then is used for the derivation of 
the meson-exchange potential in coordinate space. The formalism is developed to explicitly consider S-D wave mixing 
and also to include isospin-symmetry breaking effects. In Sec. Ill we first discuss the numerical procedure for solving 
the coupled-channel Schroedinger equation. Then we present our results for a possible binding of the DD* system, 
discussing the influence of various components both in the potential and the di-meson configuration. In Section IV 
we extend our approach to the BB* bound state system as well. Finally, in Sec. V, we present a short summary of 
our results and some general comments. 



II. METHOD 



A. Effective Lagrangian 

Our starting point (as in the potential approach pursued in Refs. [29, 33]) is the effective Lagrangian of heavy 
hadron chiral perturbation theory (HHChPT) [35, 36] based on chiral and heavy quark symmetries. It gives rise to 
the interaction Lagrangians between D(D*) mesons and light pseudoscalar, scalar and vector mesons with 

/Wp = -tgvv*Av a v;\-v; a vl)d»F ab , 

zw v = -2f VT> , Y e^ a (s(d^Y v u[CDld a vf - d a v\vf) - (V^d a V b - d a V*^V b )\ , 

Cvv a = -2m- D g a 'D a 'D'l l a , (1) 

Ct>.t>.„ = 2mv*9aK a Kl<7 > 

Cvvy = -^OTp!^-^9 M ^)(V)«t, 

C V , V , Y = igwoMK^WVlj - Z>; 6 0"Z>r + )(V M U + iifw^;] a K,b(d^ v - d v W) ab , 

where T>^ = (D°(*\ D + (*\ Df^). The octet pseudoscalar P and the nonet vector V meson matrices are defined as 

+ ^ - + K " 

«- + ^ ^ \ (2) 



and 




V= p- » K*° \ ■ (3) 



The HHChPT couplings are given as 



gvvr = -j- \J m-Dm-D* , gvvv — gvvv — 

Jtt v2 

3 = 0.59, /3 = 0.9, A = 0.56GcV _1 , / 7r = 132MeV, g K = 3.73 . 
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We assume that the X(3872) is a bound state of D*and D with quantum numbers J PC = 1 ++ . To compose a 
1 ++ state with a pseudoscalar and a vector meson, the orbital angular momentum can have the values of L = and 
2. The leading-order meson exchange diagrams constructed with the use of these Lagrangians are shown in Figs. 1 
(pseudoscalar- vector direct channels) and 2 (pseudoscalar- vector crossed channels). All diagrams are fixed to conserve 
angular momentum and parity. In the isospin symmetry limit, where and Z)(*) ± have no mass difference, the 

X(3872) can be represented explicitly as 

\X) = ^[(l^ ^* ) - c\D*°D )) ± (\D+D*-) - c\D*+D-))\ (5) 

where the ± sign is determined according to the total isospin. The + sign corresponds to the isovector state and the 
— sign to the isoscalar one; the sign c = ±1 adjusts even and odd charge conjugation parity, respectively. Here we 
use the convention where the vector field changes sign after charge conjugation, (for a detailed discussion about this 
issue and conventions used in literature see Refs. [29, 31]). 



B. Effective potentials in the direct channels 



The effective potentials in momentum space are derived for the direct channels as follows : 



1 



' q 2 — m 2 



1 



+ 



q 2 — m 2 



+ 



AmDrnD'tn 2 



q 2 - ml 



4m D m D *ml 



(6) 



where 7 = 1 and —3 corresponds to the isovector and isoscalar channel, respectively. To avoid the singular behavior at 
small distances we further regularize the potential with a form factor F(q 2 ) = ^ 2 ~™ 2 (m is the mass of the exchanged 
meson) at each vertex. The value of the parameter A is not strictly determined. It should be about 1 GeV (a typical 
scale in low-energy physics) though its specific value can depend on the particular application. A larger value of 
A enhances the potential at short distances hence possible binding energies will depend on A as we shall discuss in 
Sec. IV. 

The potentials in coordinate space can be obtained by Fourier transformation : 



V dil (r) = - Kl V a P a (r) 
with the radial dependence 



l[KiV P p {r) + K 3 V Q P (r)} - [K2V P u (r) + KaV Q u (r)] 



mi V 



A 2 



m^r 



2m 7T x 



')■ 



(7) 



(8) 



Qm{t) = 



\k 2 -m 2 M Y 1 



M 



2x' 



-XT 



(9) 



where 



Mm 



M 



{my — mi) 2 and x 2 = A 2 — (my — mi) 2 . The subscripts at P and Q in Eq. (7) refer to the 



corresponding exchanged meson of mass mu- Here we also define the constant Vq 



v _ ml g 2 



(10) 



which was already used in Refs. [11, 31]. This quantity is introduced to compare the strength of the one-pion exchange 
potential to the counterparts from the exchange of other scalar and vector mesons. The relevant strength coefficients 
are listed in the following. They are all dimcnsionless and the dimension of the potential is carried by Vq: 



Kl = 



K2 



lq 2 f 2 

5|ilL^ 0.146, 
8g z mj 



K3 



4g 2 " 

3 I3 2 



1.75, 



16 g 2 murno* 



0.070. 



(11) 
(12) 
(13) 
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C. Effective potentials in the crossed channels 

In the crossed channels we get the following effective potential in momentum space: 

Vc ( q \ = ,. c alv-F jKgj • g)(4 • g) , c gpp-p 3 ( £ i • gX4 • g) 

2Amx>mv q 2 — m 2 12mx>mx>» q 2 — m 2 

~ f 2 g B (g'l-g2')-(gl-g)(4-g) «w2 ^(gl -gy)- (fj -g)(4 -g) 

+ Ijcfvvy 2 2 2c .tvvv 2 2 V 14 ) 

q 2 -to 2 q 2 - to 2 

where e\ and e^, are the polarization vectors of ingoing and outgoing vector mesons, respectively. 

The quantity c = ±1, related to charge conjugation in Eq. (5), appears in every crossed diagram, which is different 
from the direct channels. The crossed channel potentials in coordinate space are obtained in similar manner as in the 
previous case: 

Wr) = c{ - \v [C-{r) + S 12 T"(r)] + \v Q [C+{r) + S 12 T+(r)} 

+ wV [C+(r) \s 12 T+(r)] K,V [C+(r) \s 12 T+{r)}}, (15) 

where 

T+(r) = ( A1 2 r 2 + 3 M r + 3)^ T -( X V+3 X r + 3)^ rT — (%r + l)- , (17) 

S12 = 3(ei • f)(e*, • f) - (el • ej,)- ( 18 ) 

Here to is the mass of the exchanged meson; ^ and x are defined as /^ 2 = to 2 — (m£>* — mj) 2 and x 2 = A 2 — (m£>* — mo) 2 , 
respectively; C + and T + are the potentials for positive values of [i 2 . In the case of 7r exchange, especially for 7r° 
exchange, the mass of 7r° is smaller than the mass difference {mo* — mo) so that /i 2 becomes negative. To take care 
of this case we define jl 2 = — fj 2 and obtain the real parts of these potentials as 

= _^c OS (Mr)- e -^_ X (A 2 -TO 2 ) e _ xr; 

r W = -M r +3 3-^-4-3^ 3-5-- [X r +3%r + 3) 3-3 — (xr + 1)- . (20) 

A is sufficiently large so that the sign of x 2 is not affected by the sign of fj 2 . 

The functions C and are dimensionless and their shape depends on the value of A and the mass to of the 
exchanged mesons. To compare the strength of the potential due to vector meson exchange to the counterparts of the 
pseudoscalar mesons we define the dimensionless parameter K4 as 



2 m 2 



^ = \fl*vv^t = 2A2 ^T - lm - ( 21 ) 
3 g 2 g 2 

We present the potentials in the crossed channels again to see the tensorial terms clearly. Writing the potential in 
the L = 0, 2 basis we have 



Across (r) = V 

where 



1 \ „ , . , / -V2 



Ql) Across (r) + I 1 ) Across (»") 



(22) 



Ccro-(r) = -\c-{r) + l -C+{r) + 7 « 4 C p +(r) - Ki C+(r), (23) 
TcrossM = -\T-{r) + \T+{r) - ^K 4 T+(r) + \^T+{r) . (24) 

Finally we get the total meson exchange potential, which is the sum of the potentials from the direct and crossed 
channels: 

Kotai(r) = V dir (r) + V cmss (r). (25) 
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The profile functions P(r), Q(r), C'(r) and T(r) encoding the contribution of meson exchange to the potential are 
plotted in Fig. 3. One can see that the functions C(r) and Q(r) dominate at small distances. Of course strength and 
range of all curves depend on the mass of the exchanged meson and the value of A, but Fig. 3 contains at least the 
characteristics of each function. Vb is fixed at 1.3 MeV (or 1.5 McV) by the experimental data on the decay width 
T(D* + — ► D°n + ). We choose the value Vb=1.3 MeV used before in Refs. [11, 31] which is also consistent with the 
parameters of the HHChPT Lagrangian (4). 

Multiplying the relevant coefficients such as KiVo, K2V0, etc. with P(r), Q(r), C(r) and T(r) in Fig. 4 we indicate 
the potentials in the isoscalar limit (that is, with 7 = —3). The total potential is written as 

«--%(i?)+w(_5r^) (26) 

in the L = 0, 2 basis. The tensor potential Vr — VoT CIOSS arises only from the crossed channels. The central potential 
Vc = Vdi r + VoC cr0 ss does not mix S and D waves. In Fig. 4(a) we indicate the individual contributions to Vc and 
the total result for a cutoff in the form factors of A = 1250 MeV. The contributions of n and p meson exchange are 
dominant and attractive. The effect of cr-exchange is small, though it depends on the a mass. The vector meson 
potentials become slightly attractive around 0.1 - 0.2 fm but are negligible compared to the total potential. In Fig. 4(b) 
we indicate the off-diagonal potential Vr including all exchanged mesons for A=1250 MeV. The pseudoscalar and 
vector meson potentials have different sign, but the tt contribution dominates such that the total potential becomes 
attractive. ^From Eq. (14) it is evident that the different signs of the pseudoscalar and vector meson potentials are 
dictated by the opposite sign in the effective Lagrangian. Not only Vc plays a possible role to generate attraction for 
the X state but also Vr though its strength is small compared to Vc at short distances. 



D. Isospin symmetry breaking 

We set up the A(3872) as a mixture of the isoscalar and the isovector components, that is \X) = c |0) + c±|±), 
where |0) = -±(\D°D*°) - c\D*°D )) and |±) = -j-QD+D*-) - c\D* + L>-)). In the limit of isospin symmetry the 

coefficients exactly fulfill the relation |co| = |c±| = -^=. The mass differences between the neutral and charged channel 
induce slight changes in these coefficients. When isospin symmetry is broken, the isoscalar factor 7 = — 3 in the 
potentials should be replaced by a 2 x 2 matrix 



-1 -2 
-2 -1 



(27) 



in the particle basis of neutral |0) and charged |±) states. Here the diagonal matrix elements of —1 come from the 
exchange of neutral mesons, while the off-diagonal matrix elements —2 arise from charged meson exchange. We have 
four different potentials according to the corresponding diagrams. Then the Schrodinger equation becomes as follows: 



Mi 



v 2 

2m 









-V a -2V C 
-2V d -V b 



co|0) 
c±|±) 



= E 



co|0) 
c±|±> 



(28) 



where mo and m± are the reduced masses while Mq, M± are the total masses of neutral and charged systems, 
respectively. Here, for example, V a is the potential which is generated by the diagrams (a) in Figs.l and 2 with 



V a = 



1 
1 



V 



dir 



1 
1 



V C? 







-V2 



V2 1 



(29) 



in the L — 0, 2 basis. Accordingly, Vb, c ,d & re due to the corresponding graphs of Figs.l and 2. 



III. RESULTS 



A. Numerical method 



To get solutions to the Schrodinger equation with the derived potential we need to solve coupled second-order 
differential equations. The numerical procedure is as follows: we discretize r-space and diagonalize the potential 
at each discretized position with the boundary conditions that both S and D waves vanish at r — > 00. Then the 
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Hamiltonian becomes a finite matrix and we solve the Schrodinger equation using diagonalization routines. This 
method is well suited for obtaining the ground state wave function and eigenvalue in the bound system. We tested 
the results varying the values of the boundary position and the number of discretization points. 

To confirm our numerical calculations we also adopt another method in solving the Schrodinger equation by using 
MATSCS, which is a Matlab package implementing the CPM{P,N} methods for the numerical solution of the mul- 
tichannel Schrodinger eigenvalue problem [37]. We obtain agreement in the results for the energy eigenvalue within 
0.1 MeV between our matrix method and MATSCS. This double check can only be performed in the limit of isospin 
symmetry, because MATSCS can only solve symmetric potentials. When isospin symmetry is broken the multichannel 
potentials are slightly distorted and do not coincide any more for the neutral and the charged states. 

B. Results 

In Table I we summarize our results for the binding energy Dbi n = E — M (M is the total mass of the neutral 
component) in dependence on the form factor cutoff A. In addition we also indicate the probabilities P for having 
the neutral |0) or the charged |±) components either in S- or D-wave in the bound state wave function. The size of 
the system is characterized by the rms radius of the dominant neutral S— wave component. 

As A is growing the binding energy becomes larger. This is expected since the attractive potentials get stronger 
when A is growing. The ratio between the neutral and the charged states is very sensitive to the explicit value of A. 
The isospin-breaking effect mainly comes from the mass difference between the neutral and the charged DD* system, 
that is M± — Mo ss 8.1 MeV. The corresponding wave functions at a binding energy of —0.40 MeV are shown in Fig. 5. 
The neutral states dominate and the D-wave components are negligible. In Fig. 6 we also present the corresponding 
probabilities near |£?bin| ~ 0. When the binding energy become larger, however, the effect of the mass difference 
between neutral and charged states weakens and the isospin symmetry is almost restored. For example, at a value of 
A=2500 MeV the binding energy is —555.34 MeV, which by its absolute value is larger than 8.1 MeV. For this value 
of A the probabilities P(0 s ),P(0 D ),P(±s) and p (±d) are 46.1, 4.2, 45.6 and 4.1, respectively. 

In Table II we indicate our results in the isospin limit. Compared with the full case of Table I the isospin-breaking 
effects reduce A by about 90 MeV to obtain a bound state with the same binding energy. The isovector state of DD* 
cannot make a bound state while the isoscalar state can and mixing of isovector and isoscalar components leads to a 
weakening of the binding energy. 

To identify the relevant components both in the potential and the bound DD* configuration we also looked at 
reduced variants of this approach. For example, when just keeping pion-exchangc in the potential a bound state can 
be formed for values of A > 1700 MeV. This result is consistent with the findings in Ref. [31]. Turning on in addition 
a meson exchange slightly increases the binding energy with the effect depending on the explicit mass value m a . We 
vary the value of m a from 200 MeV to 600 MeV. When we consider 7r and p meson exchange we only get a bound state 
for values of A > 1250 MeV. Further additional meson exchange components do not introduce a significant effect. 

When we turn off the charged DD* components we assume that the X meson is no more isoscalar, but a bound 
state of neutral D°D*° components only. Including in this case pion-exchangc only leads to a minimal value of 
A = 4450 MeV for the case of binding. Neglecting in addition the D-wave effect a bound state can only be generated 
for A > 5900 MeV. These results fully agree with the ones of Ref. [29] . If we include all meson exchanges but still 
neglect the charged component a bound state is obtained at A = 2050 MeV. Further neglect of the D-wave results in 
a minimal value of A = 2300 MeV to form a bound state. 

Switching off the D-wave results in a binding energy of -0.23 MeV at A = 1250 MeV. This value of -0.23 MeV 
is smaller than the value of —6.32 MeV which is obtained in the full calculation. Thus we conclude that to form a 
bound state the condition of I = plays an important role as well as tt and p meson exchanges, but S — D wave 
mixing effect is less significant. 

IV. BB* BOUND STATES 

Heavy quark symmetry and the nonrclativistic approximation is more reliable for heavy-light systems containing a 
b instead of a c quark. ^From this point of view it is worthwhile to study if in analogy to the DD* system B and B* 
mesons can also form molecular bound states. In Table III we indicate possible BB* bound states in dependence on 
A. The mass values of to#=5279.4 MeV and ms*=5325.0 MeV are used as input. The mass differences between the 
neutral and charged mesons are not known experimentally at the moment. Therefore, isospin-breaking effects 
only arise from the mass differences of the exchanged mesons such as tt or p. However, these mass differences for the 
exchanged mesons cause very slight deformation in the potential and its effect is negligible. Presently it is the best 
option to assume isospin symmetry. 
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We calculate combinations of possible BB* bound states assumed to be an isoscalar or an isovector with charge 
parity even or odd. BB* is likely to form a bound state both with J PC = 1 ++ and l" 1 in the isoscalar limit for 
reasonable values of A. In the case of DD* no bound state with J PC = l" 1 is possible for a suitable value of A. 
Compared to the DD* case the heavier mass of the BB* states leads to a reduction of momentum and hence of the 
angular momentum terms which act repulsive, and therefore binding becomes easier. 

As evident from the first part (related to J PC = 1++) of Table III binding energies even increase for reduced values 
of A < 700 McV. As a matter of fact the parameter A was introduced to regularize the attractive delta-function in 
the effective potential. The potential becomes deeper as A increases with the form factor approaching 1 when A goes 
to infinity. However, as displayed in Fig. 7(b) the potential gets deeper even when A becomes smaller. This effect 
happens when A is smaller than the mass of the exchanged meson. This might not have any significant physical 
meaning at this low A values, but it is questionable if one can simply ignore this behavior of a deepening potential at 
low A. Despite of this BB* bound states can be formed for a reasonable values of A, around 1 GeV, which is larger 
than the mass of the exchanged mesons. 

V. CONCLUSION 

In this paper we explore the possibility to generate bound states with D and D* mesons in a potential model. 
Using the HHChPT Lagrangian we construct the effective potential including isospin symmetry breaking and also 
S-D wave mixing. In this effective potential the whole light pscudoscalar and vector mesons play the role of exchange 
particles between D and D* mesons. Because of the heavy mass of the DD* system a nonrelativistic approximation 
is meaningful. We therefore solve a coupled channel Schodinger equation with specific potentials that are given by 
the effective Lagrangian. Binding of the DD* system is obtained for reasonable values of A ~ 1.2 GeV, increase of 
this cutoff value will generate deeper binding. Comparing this result to the isospin symmetry limit the DD* bound 
state where isospin symmetry is broken needs a somewhat larger A but still within a reasonable range of values. By 
switching on and off various factors we also demonstrated that the relevant ingredients for binding are the 7 = 
condition and 7r, p meson exchanges. The S-D wave mixing effect is essentially negligible. 

We extended our method to the BB* system to see if a molecular bound state is possible. Bound states of BB* 
are formed both with J PC = 1 ++ and 1 H in the isoscalar limit. But we find that it is harder to make a BB* bound 
state in the isovector limit just as for DD*. If the mass differences between the neutral and charged B^ were known, 
then we could apply our method to see if isospin breaking effects play an important role. But presently we assume 
isospin symmetry in calculating the binding energies of the BB* -system. 

It is remarkable that there could be two different states in a molecular BB* picture: with J PC = 1 ++ and V 
(both with I = 0). For example when A = 1000 MeV the binding energy of 1 ++ is —21.13 MeV while that of 1 H is 
—0.34 MeV. Future possible detection of these states can give strong support to the molecular approach in the heavy 
meson sector. 
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FIG. 1: Feynman diagrams describing the scattering of D and D* in the pseudoscalar- vector direct channel: (a) and (b) - 
Isospin channel I — and the isotopic factors are equal to —1; (c) and (d) - Isospin channel I — 1 and the isotopic factors 
are equal to 2. The bold (solid) lines stand for vector D* (scalar D) mesons while the dashed lines represent the exchanged 
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FIG. 2: Feynman diagrams describing the scattering of D and D* in the pseudoscalar-vector crossed channel; (a) and (b) - 
Isospin channel I — and the isotopic factors are equal to —1; (c) and (d) - Isospin channel I — 1 and the isotopic factors are 
equal to 2. The bold (solid) lines stand for vector (scalar) mesons while the dashed lines represent the exchanged mesons. 
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FIG. 3: Comparison of the shapes of the dimensionless potentials for A = 1250 MeV. The mass value is chosen as the p meson 
mass, rn = 771 MeV. 
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FIG. 4: Central and tensor potentials due to meson-exchange for A=1250 MeV. 

(a) the central potentials Vc (n and p meson exchange is dominantly attractive); (b) the tensor potentials Vr due to the crossed 
channels (it and p mesons are dominant and equal in sign) . 
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FIG. 6: The corresponding probabilities near |i?bin| ~ 0. Bold lines denote the neutral states while dot-dashed lines denote the 
charged states. 
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FIG. 7: The behavior of the central potentials for various values of A. 
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TABLE I: Binding energies with varying A in the isospin breaking case. 



A [MeV] 


E hin [MeV] 


E = Mo + E bin [MeV] 


P(0s)% 


P(0d)% 


P(±s)% 


P(±d)% 


rms(Os) [fm] 


1100 


No bound state 


- 


- 


- 


- 


- 


- 


1136 


-0.10 


3871.10 


91.9 


0.3 


7.5 


0.3 


7.5 


1150 


-0.40 


3870.80 


86.0 


0.5 


13.0 


0.5 


4.8 


1160 


-0.71 


3870.49 


82.1 


0.6 


16.7 


0.7 


3.7 


1168 


-1.02 


3870.18 


79.2 


0.7 


19.3 


0.7 


3.1 


1200 


-2.65 


3868.55 


70.5 


1.0 


27.5 


1.0 


1.9 


1250 


-6.32 


3864.88 


62.6 


1.3 


34.9 


1.3 


1.3 


1300 


-11.10 


3860.10 


58.2 


1.5 


38.8 


1.5 


1.0 


1350 


-16.87 


3854.33 


55.5 


1.7 


41.1 


1.7 


0.8 


TABLE II: Binding energies with varying A in the isospin symmetry limit. 


A [MeV] 


E bin [MeV] 


E = M + E hin [MeV] 


P(0s)% 


P(0d)% 


P(±s)% 


P(±d)% 


rms(Os) [fm] 


1000 


No bound state 














1050 


-0.06 


3871.14 


49.7 


0.3 


49.7 


0.3 


5.9 


1076 


-0.41 


3870.79 


49.5 


0.5 


49.5 


0.5 


3.7 


1090 


-0.71 


3870.49 


49.4 


0.6 


49.4 


0.6 


2.9 


1100 


-0.97 


3870.23 


49.3 


0.7 


49.3 


0.7 


2.5 


1150 


-2.94 


3868.26 


49.1 


0.9 


49.1 


0.9 


1.6 


1200 


-5.95 


3865.25 


48.9 


1.1 


48.9 


1.1 


1.2 


1250 


-9.98 


3861.22 


48.7 


1.3 


48.7 


1.3 


0.9 



TABLE III: Binding energies of the BB* system. 



Isospin 




A [MeV] 


E hin [MeVl 


E = Mo + £ bin [MeV] 


P(S)% 


P(D)% 


rms(S') [fm] 






600 


-6.32 


10598.08 


97.90 


2.10 


1.24 






650 


-1.67 


10602.73 


97.65 


2.35 


1.91 






700 


-0.39 


10604.01 


97.86 


2.14 


3.41 




1++ 


750 


-0.33 


10604.07 


97.72 


2.28 


3.62 






800 


-1.09 


10603.31 


96.90 


3.10 


2.15 






850 


-3.28 


10601.12 


96.32 


3.68 


1.38 






900 


-7.33 


10597.07 


95.97 


4.03 


1.02 


1=0 




950 


No bound state 














1000 


-0.34 


10604.06 


92.70 


7.30 


3.92 






1050 


-2.01 


10602.39 


86.24 


13.76 


1.84 




1+- 


1100 


-5.88 


10598.52 


82.06 


17.94 


1.20 






1150 


-12.72 


10591.68 


79.45 


20.55 


0.89 






1200 


-23.22 


10581.18 


77.83 


22.18 


0.71 






1250 


-38.00 


10566.40 


76.80 


23.20 


0.59 






4700 


No bound state 














4750 


-0.01 


10604.39 


98.98 


1.02 


9.01 






4800 


-0.06 


10604.34 


98.40 


1.60 


7.00 




1++ 


4850 


-0.14 


10604.26 


97.70 


2.30 


5.25 






4900 


-0.26 


10604.14 


96.97 


3.03 


4.03 






4950 


-0.41 


10603.99 


96.24 


3.76 


3.22 






5000 


-0.62 


10603.78 


95.52 


4.48 


2.67 


1=1 




1950 


No bound state 














2000 


-0.04 


10604.36 


99.45 


0.55 


7.47 






2050 


-0.17 


10604.23 


99.10 


0.90 


4.75 




1+- 


2100 


-0.37 


10604.03 


98.79 


1.21 


3.28 






2150 


-0.67 


10603.73 


98.53 


1.47 


2.50 






2200 


-1.06 


10603.34 


98.30 


1.70 


2.02 






2250 


-1.55 


10602.85 


98.11 


1.89 


1.70 



